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Abstract: We analyze nested Bethe ansatz (NBA) and the corresponding finite size cor- 
rections. We find an integral equation which describes these corrections in a closed form. 
As an application we considered the conjectured Beisert-Staudacher (BS) equations with 
the Hernandez-Lopez dressing factor where the finite size corrections should reproduce 
generic one (worldsheet) loop computations around any classical superstring motion in the 
AdS^ X background with exponential precision in the large angular momentum of the 
string states. Indeed, we show that our integral equation can be interpreted as a sum 
over all physical fluctuations and thus prove the complete 1-loop consistency of the BS 
equations. In other words we demonstrate that any local conserved charge (including the 
AdS Energy) computed from the BS equations is indeed given at 1-loop by the sum of 
charges of fluctuations up to exponentially suppressed contributions. Contrary to all pre- 
vious studies of finite size corrections, which were limited to simple configurations inside 
rank 1 subsectors, our treatment is completely general. 
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1. Introduction 



Bethe equations describe the scattering of the fundamental degrees of freedom of inte- 
grable 1 + 1 dimensional theories defined on some large circle of length £. The existence of 
a large amount of conserved charges results in the factorizability property of the scattering 
matrix. Namely the full n particle 5-matrix is completely fixed by the 2 particle scat- 
tering. Moreover in two dimensions this 2 to 2 scattering process conserves not only the 
total momentum but also the set of individual momenta. Then, for a large enough circle, 
the momenta of the several particles are quantized through the wave function periodicity 
condition 

L 

l = e^P>^^l[S{pk,Pj) (1.1) 

meaning that the (trivial) phase acquired by a particle with momentum while going 
around the circle equals the free propagation plus the scattering phases shifts (or time 
delay in coordinate space) due to the passage through each of the other particles. In 
general Bethe equations are only asymptotically exact as £ — > oo otherwise wrapping 
effects |2|, ^ must be taken into account. 



Equation (1.1) is, however, describing particles with no isotopic degrees of freedom, 
that is S{pk,Pj) is just a phase. In general, when we have some nontrivial symmetry group, 
this is not the case and, rather, we must solve the diagonalization problem 

L 

|V') = e*^'=^n S{pk,Pj)\^l^) 

where S{pk,Pj) is now a matrix and is the multi-particle wave function (for integrable 
theories the number of particles is conserved). If the scattered particles transform under 
some symmetry group we will obtain not just one equation like ( |1.1| ) but rather a set of 
n equations entangling the scattering of particles with momenta pk and pj in space-time 
with the scattering of spin waves in the isotopic space. 

In this paper we will mainly consider the particular limit of low energies when the 
wave length of the spin waves are large and particles exhibit collective behavior which, in 
some important cases, can be associated with the classical motion of collective fields. By 
studying carefully this limit one can get important information about the quantization of 
some classical field theories. 

In terms of the Bethe ansatz equation this corresponds to a limit, first considered in 
the condensed matter literature by Sutherland |Q] in the study of the ferromagnetic limit 
of the Heisenberg chain and rediscovered and generalized in the context of AdS/CFT Q, 
where the Bethe roots Uj ~ cot(pj/2) scale with the number of such roots and with the total 
number of particles, uj ~ Ka ~ L. In this limit the Bethe roots condense into disjoint cuts. 
Since there are several types of Bethe roots, one for each Bethe equation, the condensation 
of the Bethe roots for systems with n Bethe equations will generate some Riemann surface 
with n + 1 sheets as in figure 0. This resulting curve is in one-to-one correspondence with 
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the curves classifying classical solutions through the finite gap method 0, |8|, 1C|. In 
this way one finds the semi-classical spectrum of the theory. 

The next natural step is to compute the first quantum corrections to the semi-classical 
spectrum, which from the Bethe ansatz point of view will correspond to the finite size (i.e. 
1/L) corrections. For the simplest Bethe equations of the form (1.1) these corrections, 
called anomalies, were known [0, |l2|, |l|, 0, ^, |l|, ^ but for nested Bethe ansatz 
equations the analysis is much more delicate due to the formation of bound states, called 
stacks ||l^ , which are the basic constituents of the cuts made out of more than one type of 
Bethe roots like the ones in figure ^. In this paper we develop the necessary tools to deal 
with these richer systems with isotopic degrees of freedom. 

Particularly important tools are the so called dualities. One of them, the fermionic 
duality, is well studied |19, 2C, 21, 22, 24 1 and has a clear mathematical meaning. If 

the symmetry group under which the fundamental particles transform is a super group then 
there are several possible choices of NBA equations corresponding to the several possible 
choices of super Dynkin diagram which, for super-groups, is not unique. These equations 
are related by some dualities associated with the fermionic nodes of the corresponding super 
Dynkin diagram. In the scaling limit they correspond to the exchange of Riemann sheets. 
In this paper we also use an analogue of this duality, baptized bosonic duality, which exists 
even in the case of a purely bosonic symmetry. It is associated with the bosonic nodes of 
the Dynkin diagram. 

Next we apply our method to the recently conjectured Beisert-Staudacher (BS) Bethe 
equations ||2^ . These equations contain a free parameter A and should describe two systems 
at the same time: four dimensional = 4 SYM and type IIB super-strings in AdS^ x S^, 



two theories which are conjectured to be dual ||2^, 27, 28]. At weak coupling, A <C 1, we are 
in the perturbative regime of = 4 SYM and the Bethe equations describe the spectrum 
of the planar dilatation operator which can be considered as a spin chain Hamiltonian 
|29, 30| with PSU{2,2\A) symmetry. At strong coupling \/A L ^ 1 the theory describes 
classical super-strings in the curved space- time AdS5 x |, 0] and the 1/VX corrections 
in the scaling limit correspond to the semi-classical quantization of such highly non-trivial 
quantum field theory. 

As it was stressed in our previous papers 32 1 there are two completely different 
ways to compute the 1-loop correction to the quasi-classically quantized energy spectrum. 
One, straightforward but technically more involved, is to take the NBA equations, compute 
its spectrum and then expand it in powers of i.e. find its finite size corrections. 

Another way, more indirect one, is to pick some classical solution satisfying the semi- 
classical quantization condition, and quantize around it, i.e. find the spectrum of all 
possible excitations of this solution. The one loop shift is then given by the zero energy 
oscillations and is equal to half of the graded sum of all excitation energies, like for a simple 
set of independent one dimensional harmonic oscillators. 

Both calculations can be performed using the BS equations and it is a very nontrivial 
test of the proposed equations that these two calculations give the same result. In fact for 
the second calculation we do not even need the Bethe ansatz, since it is based completely on 
the semi-classical quantization which, as shown in [31|, can be performed relying uniquely 
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on the classical integrable structure of the theory - the algebraic curve ^ . Moreover we 
expect the second approach to give the exact result whereas the first one is only valid as long 
as one can trust the asymptotic BAE, which suffers from the wrapping effects mentioned 
above. Indeed we found that the two results coincide not precisely, but only for large 
L/^fX with exponential precision. This is obviously a manifestation of the wrapping effects 



considered in the AdS/CFT context in 
first observed in [031 . 
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Figure 1: For su{K\M) super algebras the 
Dynkin diagram is not unique. The several possi- 
ble choices can be represented as the paths going 
from the up right corner (M, K) to the origin al- 
ways approaching this point with each step. The 
turns are the fermionic nodes whereas the straight 
lines correspond to the usual bosonic nodes. Dif- 
ferent paths will correspond to different sets of 
Bethe equations which are related by fermionic 
dualities which flip a left-down fermionic turn 
into down-left turn or vice- versa p4|. 



35, 36, 37|. This exponential mismatch was 

Finally we should stress that we follow 
a constructive approach. That is we start 
from the classical integrable structure, the 
finite gap curves. The curves can be de- 
scribed by some integral equations. We find 
how to correct this equations in such a way 
that they will now describe not only the 
classical limit \/A — > oo but also the l/^f\ 
corrections. Then we show that the inte- 
gral equations modified in this way coincide 
precisely with the scaling limit expansion of 
the BS equations ||2^ with the Hernandez- 
Lopez phase |3^ (up to some exponentially 
suppressed wrapping effects, irrelevant for 
large angular momentum string states)! Our 
comparison, being done at the functional 
level, is completely general. 

This paper is organized as follows. In 
section ^ we introduce some notations, the 
notion of stack and the bosonic duality. In 
section]^ we derive, in two independent ways, 
an integral equation describing the finite 
size corrections to the leading limit - using 
the dualities and using the transfer matri- 
ces. In section^ we follow the constructive approach mentioned above to re-derive the same 
integral equation from the equations in the scaling limit. Section ^ contains some details 
about the bosonic duality such as some theorems and examples - the reader interested only 
in the main results of the paper can skip this section. In section ^ we apply the methods 
developed in the previous sections to the study of the BS equations, compute the finite 
size corrections and relate them with the quantum fluctuations of the theory. Appendix A 
is devoted to the study of the invariance of the transfer matrices of su{K\M) supergroups 
under the bosonic duality and in Appendix B we derive an integral equation describing the 
semi-classically corrected equations for su{n) spin chains. 
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2. Nested Bethe Ansatz and Bosonic Duality 



In the first sections we stick mainly to the simple example of su{l, 2) spin chain, although 
our main motivation comes from its application to AdS^ x /J\f = 4: SYM correspondence 
where the symmetry group is PSU{2,2\4:). Indeed this simple toy model contains already 
all the nontrivial new features appearing due to the Nested nature of the Bethe ansatz. The 
generalization to other (super) groups is straightforward and in particular we shall focus on 
the Bethe ansatz describing the superstring in AdS^ x in section ^. 

For integrable rank r spin chains each quantum state is parameterized by a set {ua,j} 
of Bethe roots where a = 1, . . . , r refers to the Dynkin node and j = 1, . . . , Ka where Ka 
is the excitation number of magnons of type a. The Bethe equations from which we find 
the position of these roots are then given by 



Ua,j 



Wa J Qb [Ua,j - ^Mah) 



where 



Qa{u) = W{u- Ua,j) 

i=i 

are the Baxter polynomials, Va are the Dynkin labels of the representation considered and 
Mab the Cartan matrix. In fact, contrary to what happens for the usual Lie algebras, for 
super algebras the Dynkin diagram (and thus the Cartan matrix) is not a unique. Take 
for example the su{K\M) super algebra. The different possible Dynkin diagrams can be 
identified |^] as the different paths starting from (M, K) and finishing at (0, 0) (always 
approaching this point with each step) in a rectangular lattice of size M x A' as in figure 
|l[ The turns in this path represent the fermionic nodes whereas the bosonic nodes are 
those which are crossed by a straight line - see figure |^ (the index a goes along the path 
as indicated). The Cartan matrix Mab is then given by 

Mab = {Pa + Pa+l) ^ab - Pa+l5a+l,b " Pa^a,b+l 

where pa is associated with the link between the node a and a + 1 and is equal to +1 (—1) 
if this link is vertical (horizontal). 

Here we are considering twisted (quasi-periodic) boundary conditions which, from an 
algebraic Bethe ansatz point of view corresponds to the diagonalization of a transfer matrix 
with the insertion, inside the trace, of an additional diagonal matrix [^] which can be 
parameterized by 

5 = diag {e''f'\...,e"t"^,e''P\...,e''pA e SU{K\M) (2.2) 



and the twists Tq, appearing in (2.1) and associated to a Dynkin node located at (m, k) in 



the M X K network depicted in figure |l], are then given by [40| 

Ta = 4'k ~ (pk+i for a bosonic along a vertical segment of the path 

Ta = Vm+i — for a bosonic along a horizontal segment of the path 

Ta = Vm+i — <Pk + for a fermionic node in aT like turn that is with Pa~i = —Pa = 1 

Ta = 4>k+i — + for a fermionic node with Pa-i = —Pa = — 1 
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2a 



2b 



2c 



Figure 2: The middle node Bethe roots U2 can condense into a line as depicted in figure |p (The 
spins in this spin chain transform in a non-compact representation and tlius the cuts are topicaUy 
real. For the su(2) Heisenberg magnet the solutions are distributed in the complex plane as some 
umbrella shaped curves [||.). Roots of different types can form bound states, called stacks iQ, as 
shown in figure ||b. The stacks behave as fundamental excitations and can also form cuts of stacks 
as represented in figure 



Notice that since g G SU{K\M) we have Ylk^k — Ylm^^ ~ ^ mod 27r. We shall study 
these Bethe equations with generic twists and we will see that the usual case (T(j = 0) is in 
fact quite degenerate. 

As mentioned above, we find already all the ingredients we will need for the study 
of the BS equations in the simple example of a sm(1,2) spin chain in the fundamental 
representation described by the following system of NBA equations^ 



Qi {uij - i) Q2 (uij + i/2) 



j = l...Ki (2.3) 



\U2,j + §/ Q2 [U2,j - t) Ql {U2,j + V2) 

The eigenvalues of the local conserved charges are functions of the roots U2j only and are 
given by 

i / \ 1 \ 

We will often denote these roots carrying charges by middle node roots. 

First, consider only middle node excitations, Ki = 7^ K2 in the Sutherland scaling 
limit Q where u K2 '-^ L ^ 1. We shall always use Xaj = Uaj/L to denote the rescaled 
Bethe roots in the scaling limit. Then, the Bethe equations in log form, to the leading 
order, can be cast as 

^2,j L ^ X2,j - X2^k 

where the integers rij come from the choice of the branch of the logs. 



^These equations are exactly the same as for the su{3) spin chain except for the sign of the Dynkin 
labels which makes the system simpler because the Bethe roots are in general real. 
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3a 3b 



Figure 3: In the scaling limit, to the leading order, the bosonic duality reads Q2 — QiQi with 
Qa — Y\^:^i{u ~ "a)- Thus, if we start with the configuration in figure ^ where the Ki roots ui 
form a cut of stacks together with Ki out of the K2 middle node roots U2 and apply the bosonic 
duality to this configuration, the K2 — Ki new roots ui must be close to the roots U2 which were 
previously single while the cut of stacks in the left of figure ^ will become, after the duality, a cut 
of simple roots - see figure ^d. 



We see that we can think of the Bethe roots as positions of 2d Coulomb charges on a 
plane with an external potential equal for every particle plus an external force 27rnj specific 
of each Bethe root. Thus, if we group the K2 mode numbers {%•} into large groups of 
identical integers and consider the limit where both L and K2 are very large with K2/L 
fixed, the Bethe roots will be distributed along (real) cuts C^, each parameterized by a 
specific mode number {n"^} where A= 1, . . . ,N. Then the equations can be written 
through the density of middle node roots X2 as 

27rn^ + (p2-<p3 = - + 2(p2ix), x G (2.7) 

X 

where we introduce the resolvents 

Ga{x) = [ ^ , pa{y) = i V 5(X - Xa,j) (2.8) 

X — V L ^-^ 

and where the slash of some function means the average of the function above and below 
the cut, (J?(x) = ^ {G{x + ie) + G{x — ie)). Let us also introduce some notation useful for 
what will follow. Defining the quasi-momenta as 

Pi = -— +Gi -01, 

P2 = -^-Gi+G2-02, (2.9) 

P3 = -77- - ^2 - 03 , 

2x 



we can add the indices 23 to the mode number and to the cut C"^ in (2.7) and recast 
this equation as 



27rn^3 = i/2 - , a; G C2^3 . (2.10) 
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Next let us consider a state with only two roots U2.1 = u and ui^i = v with different 
ffavors, that is Ki = K2 = 1. Bethe equations then yield 



1 .4'i-(t>-i + 27rn I 4>i-(t)2 , . 

u = - cot , V = u-\ — cot (^-ll) 

which tells us that if n ~ 1 we are in the scaling limit where f ~ li ~ L and v = u + 0{1) 
- the two Bethe roots form a bound state, called stack and can be thought of as a 
fundamental excitation - see figure 2b. On the other hand we notice that, strictly speaking, 
for the usual untwisted Bethe ansatz with (j)a = the stack no longer exists. 

Since the stack in figure 26 seems to behave as a fundamental excitation one might 
wonder whether there exists a cut with Ki = K2 roots of type ui and U2, like in figure 2c, 
dual to the configuration plotted in figure 2a. To answer affirmatively to this question let 
us introduce a novel kind of duality in Bethe ansatz which we shall call bosonic duality. 

Indeed, as we explain in detail in section |5|, given a configuration of Ki roots of type 
ui and K2 roots of type U2, we can write 

2i sin (r/2) Q2iu) = e"l'^Qx{u - i/2)Qi{u + i/2) - e'^/^Qiiu + i/2)Qi{u - i/2) , (2.12) 

where 

Ki 

Qi{u) = n - ^ij) , Ki = K2-Ki, 

and T = (pi — (p2- Moreover this decomposition is unique and thus defines unambiguously 
the position of the new set of roots ui. Then, as we explain in section ||, the new set of 



roots {ui,U2} is a solution of the same set of Bethe equations (2T) with 



^2 • 



Let us then apply this duality to a configuration like the one in figure 2a where the roots 
U2 ^ L are in the scaling limit and where there are no roots of type ui, Ki = 0. To the 
leading order, we see that the ui in (2.12) will scale like L so that the ibz/2 inside the 



Baxter polynomials can be dropped and we find Q2 — Qi-, that is 

= U2,j + 0(1) 

and therefore we will indeed obtain a configuration like the one depicted in figure 2c. 
Moreover the local charges (|2.5| ) of this dual cut are exactly the same as those of the 
original cut 2a since they are carried by the middle node roots U2 which are untouched 
during the duality transformation. 

Finally, if we apply the duality transformation to some configuration like that in figure 
3a in the scaling limit we find, by the same reasons as above, that Q2{u) — Qi{u)Qi{u). 
This means that the dual roots ui will be close to the roots U2 which are not yet part of a 
stack - the ones making the cut in the right in figure 3a. Thus, after the duality, we will 
obtain a configuration like the one in figure 36. 
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4a 4b 

Figure 4: In the scaling limit the configurations in figure ^ condense into some disjoint segments, 
cuts, and we obtain a Riemann surface whose sheets are the quasi-momenta. In this continuous 
limit the duality corresponds to the exchange of the Riemann sheets. 

We conclude that, in the scaling limit with a large number of roots, the distributions of 
Bethe roots condense into cuts in such a way that the quasi-momenta pi introduced above 
become the three sheets of a Riemann surface, see figure 4a, obeying 

27rnjJ = - , a; G Cif . (2.13) 

when X belongs to a cut joining sheets i and j with mode number . The duality trans- 
formation amount to a reshuffling of sheets 1 and 2 of this Riemann surface^ so that a 
surface like the one plotted in figure 4a transforms into the one indicated in figure 46. 



3. Anomalies — finite size correction to Nested Bethe Ansatz equations 



In this section we will study the leading 1/L corrections to the scaling equations ( 2.13| ). 
Moreover since the charges of the solutions are expressed through middle node roots U2 and 
since these roots are duality invariant it is useful to write the Bethe equations in terms of 
these roots only. Let us then consider a given configuration of roots condensed into some 
simple cuts C23 and some cuts of stacks C13. Then, to leading order, at cuts C23 we have 



1 

- + 2 

X 



f P2{y)dy ^ r P2{y)dy 
1 x-y J 



27rn4 + 



23 



(3.1) 



C23 



Cia 



because in a cut C13 we have pi ~ p2 + C (1/-^^)- To study finite size corrections to this 
equation two contributions must be considered. On the one hand when expanding the self 
interaction we get 12, |l^, 14, 15, 16] 



^ilog 



U2,k - U2,j - % 
U2,k - U2,j + i 



2/ + 2 [ 

J x-y J 

C23 Cia 



P2{y)dy 1 , 

1- —irpo cot 1Tp2 

X — y L 



^As we shall see in the next section this interpretation can be made exact, and not only valid in the 
scaling limit. 
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where the 1/L correction comes from the contribution to the sum from the roots separated 
by 0{1). On the other hand the auxihary roots appear as^ 



^ilog 



U2,k - uij + i/2 



U2,k - ui,j - i/2 J X - y 

^ ' Cl3 



pi{y) 



dy = - [ ^dy- I PlM^PMdy 
J x-y J x-y 

Cl3 Cl3 



where the last term accounts for the mismatch in densities in cuts C13 and is clearly also a 
0{1/L) effect. Bellow we will compute this mismatch and find 



Pi{x)- P2{X) 



A coti 2 cot 21 — cot 



"23 



2mL 



2mL 



x€C 



13 



where A/ = f{x + iO) — f{x — iO) and 



cot - = ^1::^ cot 
^2 2 



(3.2) 



(3.3) 



Thus we find, for x £ C23, 



1 

- + 2 

X 



2 f P2{y)dy ^ r P2{y)dy 

J x-y J 

C23 C13 



y 



2^n4 



-F 



cot 



23 



Cl3 



A cot 12 dy 
x-y 2711 



(3.4) 



As explained before, if we apply the duality transformation, cuts C23 become cuts C13 and 
vice- versa and, to leading order, pi <-> P2- Thus for cuts C13 we find precisely the same 
equation (|3.4|) with 1 •i-^ 2, so that for x G C13 



1_^2-/" ^^M^_^ / P2{y)dy 
X J x-y J 

C13 C23 



27rn^3 + ^ 



1 

I 



cot 13 



C23 



A cot 12 dy 
y 27ri 



X 



(3.5) 



These two equations describing the finite size corrections for the two types of cuts of the 
stt(l,2) spin chain are the main results of this section. 

In what follows we will derive this result from two different angles. Namely, we will 
find this finite size corrections using a Baxter like formalism based on transfer matrices for 
this spin chain in several representations and by exploiting the duality we mentioned in 
the previous section. It will become clear that the generalization to other NBA equations 
based on higher rank symmetry groups is straightforward. 

3.1 Derivation using the transfer matrices 

The central object in the study of integrable systems is the transfer matrix T{u). The 
algebraic Bethe ansatz formalism has the diagonalization of such objects as main goal and 
the Bethe equations appear in the process of diagonalization (see |41] and references therein 
for an introduction to the algebraic Bethe ansatz). As functions of a spectral parameter u 
and of the Bethe roots Uaj these transfer matrices seem to have some poles at the positions 



^recall that the Bethe roots U2,k belongs to a C23 cut and therefore is always well separated from uij 
roots which always belong to C13 cuts. 
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of the Bethe roots. On the other hand they are defined as a product of R operators which 
do not have these singularities. This means that the residues of these apparent poles must 
vanish. These analyticity conditions (on the Bethe roots) turn out to be precisely the Bethe 
equations, and thus, if we manage to obtain the eigenvalues of the transfer matrices, we can 
use this condition of pole cancellation to obtain the Bethe equations without going through 
the algebraic Bethe ansatz procedure, see for example |42, For the su{l, 2) spin 



chain we have the following transfer matrices in the anti-symmetric representations: 

(3.6) 



Qiiu+i) Q2{u-i' 




Qi(n + i) \u-^) g2(n-i 

Tg (u) = fa (u) 

One can easily see that the Bethe equations do follow from requiring analyticity of these 
transfer matrices. 

In [^] it was shown and emphasized that the TQ relations are the most powerful 
method to extract finite size corrections to the scaling limit of Bethe equations. 

In this section we will use the transfer matrices presented above along with the fact 
that, due to the Bethe equations, they are good analytical functions of u to find what are 
the finite size corrections to this Nested Bethe ansatz. Since for generic (super) nested 
Bethe ansatz the transfer matrices in the several representations are known, this procedure 
can be easily generalized for other NBA's. 

The key idea to find the finite size corrections to NBA is to use the transfer matrices in 
the various representations to define a new set of quasi-momenta Qi as the solutions of an 
algebraic equation whose coefficients are these transfer matrices. For example, to leading 
order, 

Tn (n) ~ e'P^ + e'P^ + e'P^ , 

(^U'j ~ g«(Pl+P2) _|_ gi(p2+P3) _|_ gi(P3+Pl) ^ 

Tg (n) ~ e*(f^+P^+P3) , 
so that if we define a set of exact quasimomenta qi hy'^ 

Tg in) - e^^ Tg (n) (l " ^) + e'^^ Tu (u) (l " ^) " e''' = , (3.7) 





^Exploiting the similarity between this definition equation and 4.1 in |24| we can easily generalize this 
algebraic equation to a more general su{K\M) super group. More precisely we identify e^^" <-» e'' which 
is obviously natural if we look at 4.2 in this same paper (see also Appendix A where we use this two 
expressions slightly modified to match our normalizations). We thanks V.Kazakov for pointing this out to 
us. 
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then, to leading order, — Pi- Notice however that the coefficients in this equation have 
no singularities except some fixed poles close to u = 0. Thus, defined in this way, the 
quasi-momenta qi constitute a 4 sheet algebraic surface (modulo In ambiguities) such that 

(^i - = 2Trnfj , X e Cij , (3.8) 

and, needless to say, this is an exact result in L, it is not a classical (scaling limit) leading 
result like (2.13)! On the other hand, the expansion at large L of the above algebraic 
equation yields 

qi=Pl + ^(+COti2 +COti3) 

1 / 

q2=P2 + ^(-COt21 +COt23) 
1 / 

q3=P3 + ^(- COtsi - COt32) , 



which follows from the expansion 
To (n) (l - = e^Pi + e'P^ + e'P^ 



^ [e^'P- {2p[ - p', - p',) + e^f^ {p[ - p',) + e'P^ ip[ + p', - 2p',)] + O (-1) 



4L 



1 

"iZ 



Tg (n) = e*(Pi+P2+P3) + O (^-^ 

of the several transfer matrices. Then, to the first order in 1/L the exact equation ( p.^ ) 
gives, for the quasi-momenta pi introduced in (|2.9| ), 

i2-iz = 27m^3, - cot23 , a; G C23 (3.9) 

1^1 - = 27rn^3 - ^ (coti2 +2 cot 13 + COt32) , X eCi3 (3.10) 

where in ( |3.9| ) we use the fact that function cot3i — cot2i vanishes under the slash on the 
cut C23 since 

cot+. = cot^^. , X GCik- (3-11) 

Equations ( |3.9D ,( 3?TC| ) are the finite size corrections we aimed at! 
Finally q2 must have no discontinuity at a cut C13 and therefore 

Ap2 = 27rz (pi - /)2) = -^(cot^i -cot^g) , x G C13 . (3.12) 



Thus, replacing the quasi-momenta pi by its expressions in terms of resolvents ( |2.9D and 
relating the density of auxiliary roots pi to that of the middle node roots p2 through ( 3.12| ), 
we recover precisely (^) and ( |3.5[ ) as announced. 
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We would like to stress the efficiency of the TQ relations. We were able to find the usual 
cot contributions (coming from the expansion of the log's of the Bethe equations when the 
Bethe roots are close to each other) plus the mismatch in densities of the different types of 
roots making the cuts of stacks using only the fact that due to Bethe equations the transfer 
matrices in several representations were analytical functions of u. The computation done 
in this way is by far more economical than a brute force expansion of the Bethe equations. 

Finally let us make an important remark. To derive ( |3.5D from ( p. 10 ) one should use 



1 f ^ coti2 
coti2 = -— / dy (3.13) 

Zm J X — y 

C13UC23 

which is clearly a valid relation if coti2 has only branch cuts as singularities. For generic 
twists and for small enough cuts C13 and C23 this is the case. Indeed, in the absence of 
Bethe roots we have no cuts at all and thus pi — P2 = (^2 — (f'l- Suppose (t)2 — (pi 27rn. 
Then, by continuity, when we slowly open some cuts C23 and C13 then pi — p2 will start 
taking positive values around 02 — without ever being zero. Thus, if the cuts are small 
enough we will never get poles in cot 12. In the next section we will see that the stacks 
as described in only exist when this assumption of absence of poles is right and are 
destroyed when pi — p2 reaches 27rn. 

3.2 Re-derivation using the bosonic duality in the scaling limit 



In this section let us re-derive the mismatch formula (|3.2| ) using the bosonic duality (^.ij). 
Besides the obvious advantage for what concerns our comprehension of having a second 
derivation there are systems for which the Bethe equations are known but the algebraic 
formalism behind these equations is still not well developed (this is the case for example for 
the AdS/CFT Bethe equations proposed by Beisert and Staudacher which we will study 
in section |^ . 
Denoting 

uii = - ei , = U2A - e j , e ~ 1 



and expanding the bosonic duality (5.1) in the scaling limit (L 00) we get 



sin(r/2) = sm (i (^i - + r)) exp ( ^ + ^ 



1=1 ' 1=1 ' 

where r = 0i — 02- Taking the logarithm of this equation and differentiating with respect 
to u we get 

Ee, ii _ G'l — G'l Gi — Gi + r 

{u-ujf ^ {u - ujf ~ 2L 2 

where we notice that the left hand side is precisely the difference of resolvents G2 — Gi — Gi ! 
Thus we find 

^ ^ ~ G; - G[ Gi-Gi+T G^2 - 2G; G2-2Gi+r 1 
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Finally, by computing the discontinuity of this expression at the cuts C13 we will get the 
mismatch of the densities of the roots in a cut of stacks^ 

A COti2 cot^j^ — cot^3 
27riL 27riL ' 

which was the gap in the chain of arguments presented in the beginning of the section ^ 
and leading to (^). 

Finally let us show that the bosonic duality amounts to a simple exchange of Riemann 
sheets in the scaling limit. Consider for example 

Pi = + Gi - 4>i = --^ + G2 - Gi - (f>i = p2 
since, as we will see more carefully in the next section, (pi^2 = <p2,i- 
4. 1-loop shift 



In [31 1 we explained how to obtain the spectrum of the fluctuation energies around any 
classical string solution using the algebraic curve by adding a pole to this curve. In particu- 
lar we reproduced in this way some previous results 46, where the semi-classical 



quantization around some simple circular string motions were considered by directly ex- 
panding the Metsaev-Tseytlin action around some classical solutions and quantizing 
the resulting quadratic action. Using the fact that one extra pole in the algebraic curve 
means one quantum fluctuation, we can compute the leading quantum corrections to the 
classical energy of the state from the field theory considerations using the algebraic curve 
alone, as we mentioned in the introduction. This implies a nontrivial relation between fluc- 
tuations on algebraic curve and flnite size corrections in Bethe ansatz as we will explain in 
greater detail below. In this section we study this relation on the example of the su{l,2) 
spin chain and then in section 6 we extend this to the super-string case. 

As mentioned in the introduction, in the scaling limit u ^ K ^ L Iwe are describing 
some slow and low energetic spin waves, 

K K 

^ ' ^ U2,j + 1/4 ' ' 
around the ferromagnetic vacuum of the theory. In this limit the theory is well described by 



a Landau-Lifshistz model which is a field theory with coupling 1/L ||5^, 51, |5^. Therefore a 
very nontrivial property relating fluctuations and flnite size corrections in this NBA should 
hold: 

• Suppose we compute the energy shift 5£n due to the addition of a stack with mode 
number n uniting sheets pi and pj to a given configuration with some finite cuts C. 

• Suppose on the other hand that we compute 1/L energy expansion 
... of the configuration with the finite cuts C. 



^A/ = /+-/-, so that p= 
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From the field theory point of view the first quantity corresponds to one of the fluctuation 
energies around a classical solution parameterized by the configuration with the cuts C 
whereas the second quantity, £^^\ is the 1-loop shift around this classical solution with 
energy S^'^h This 1-loop shift, or ground state energy, is given by the sum of halves of the 
fluctuation energies |Q 

^^'^ = ^EE^c' (4-1) 

n ij 

In fact for usual (non super-symmetric) field theories this sum is divergent and needs to be 
regularized. We will see that (|4.1|) can be generalized and holds for arbitrary local charges 



2^'^ = ^EE^2^- (4.2) 



2 

n ij 

Let us stress once more that from the Bethe ansatz point of view these quantities are 
computed independently and there is a priori no obvious reason why such relation between 
fluctuations and finite size corrections should hold. In this section we will show that 
Nested Bethe Ansatz's describing (super) spin chains with arbitrary rank do indeed obey 
such property with some particular regularization procedure (for the Heisenberg su{2) spin 
chain a similar treatment was carried in |l^). Moreover we will see that the regularization 
mentioned above also appears naturally from the Bethe ansatz point of view as some 
integrals around the origin. 

4.1 1-Ioop shift and fluctuations 

In this section we will understand the interplay between fluctuations and finite size cor- 
rections in NBA's in the scaling limit. For simplicity we are considering the su(l,2) spin 
chain described in the previous sections. General su{N) is considered in Appendix B. 

Let us pick the leading order integral equation for the densities of the Bethe roots in 
the scaling limit ( |3.1| ) and perturb it by a single stack, connecting with pi. According 
to (2^) this means simply implies P2 ^ P2 + T^i^ ~ x^^)-, where x'j is position of the new 



stack. Finally, the positions where one can put an extra stack, as it follows from the BAE 



(2.3,2.4), can be parametrized by one integer mod number n 

Pi(xjj)-Pj(xjj) = 27rn. (4.3) 
Therefore, for i = 2, j = 3 the perturbed equation (|3.1D reads 



i + 2/i^+ / J!^+'^ = 2nk,, + <P,-^, , XEC23. (4.4) 
X J X — y J X — y L X — xff 

C23 C13 

and this perturbation will lead to some perturbation of the density Sp{y), which will lead 
to the perturbation in the local charges ( |2.5| ) as 



^C-^/M^>* + ^, ,4,5) 

the local charges of the fluctuation with polarization 23 and mode number n. 
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Thus, by linearity, if we want to obtain the 1-loop shift (^^) (or rather a large N 
regularized version of this quantity where the sum over n goes from — to A^) we have to 
solve the following integral equation for densities 



1 + 2/^+/ 

X J x-y J 

C23 C13 



pjy) I _}_ 

x-y ^ 2L 

^ n=-N 



X — X 



23 



+ 



1 



and then the 1-loop shifted charges are given 



Qt 



N 



^ 2L 

n=-N 
N 



C13UC23 



1 



X 



1 



13 



2'Kk 



23 



23 ■ 



(4.6) 



13^r 



ir, y n=-N V ''^ 



C13UC23 



/'^23\r ' 

cot23 <iy 



+ 



cot 13 dy 
13 y^ 2iTi 



(4.7) 
(4.8) 



To pass from the first line to the second in the above expression we use that coty has poles 
at Xn with unit residue. We will now understand how to redefine the density in such a 
way that the second term is absorbed into the first one. We start by opening the contours 
in ( [4.g| ) around the excitation points x^- These contours will then end up around the 
cuts Cki of the classical solution and around the origin. We will not consider the contour 
around x = Q - this contribution would lead to a regularization of the divergent sum in 
r.h.s. of ( ^.21 ). We will analyze it carefully in the super-string case, where it leads to the 
Hernandez-Lopez phase factor. Then we get 



C13UC23 



y'- ^ 2L 



cot23 dy 
Ci3 y'' 2vrf 



+ 



cot 13 dy 



C23 y 



2-ni 



Noting that 



cot,^. 



cot 



kj 



X £ C 



ik 



(4.9) 



(4.10) 



where the superscript -|- (— ) indicates that x is slightly above (below) the cut, we can write 



Qr 



C13UC23 



p{y) 
y' 



dy 



1 

2L 



C13UC23 



A cot 12 dy 
y"^ 27ri 



(4.11) 



so that we see that it is natural to introduce a new density, "dressed" by the virtual 
particles, 

(4.12) 



12 



^ 2L 2m 

so that the expression for the local charges takes the standard form 



Qr 



C13UC23 



y'' 



dy. 



Let us now rewrite our original integral equation (|4.6| ) in terms of this dressed density. 
We will see that the integral equation we are constructing for this density by requiring a 
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coti3 dy 
'c23 X — y 2TTi 



/■ ZXco 

J X — 

C23 



A coti3 dy 
y 2711 



C^ti3 



G 







Figure 5: Illustration of an identity used in the main text. 



proper semi-classical quantization will be precisely the equation ( p.4D which is the finite 
size corrected integral equation arising from the NBA for the spin chain! This will thus 
prove the announced property relating finite size corrections and 1-loop shift. Consider for 
example the first summand in ( [4. 61) (recall that x G C23), 



E — 



X — X. 



23 



E 



cot23 dy 
23 X — y 2'Ki 



cot 



23 



cot 23 dy 
Ci3 x-y2TTi 



cot 



23 



A cot 12 dy 

'ci3 ' 



(4.13) 

Note that cot23 has branch cut singularities at C13 which we have to encircle when we blow 
up the contour, which leads to the second term. The first term comes from the pole at 
X = y. Finally, to write the second term as it is we used ( 4.10 ). Analogously (see figure ^ 
for a pictorial explanation of the second equality) 



E— 



X — X 



13 



cot 13 dy 



Cptl3 + 



/ 

C23 



A cot 13 dy 
X — y 2-Ki 



cf^^tis - -j 



C23 



A cot 12 dy 
X — y 2711 



■ (4.14) 



Then we note that (see ( |3.13|) ) 

C)6ti3 = C^ti2 



C13UC23 



A cot 12 dy 
X — y 2711 



so that (|]^) reads 



X J X — y J X 

C23 C13 



pjy) 1 
y 2L 



2 cot 



23 



2 , Acoti2 dy _ 2 /" Acoti2 dy 



C23 



X — y 27Ti 



Cl3 



X — y 27Ti 



27r/c23 + 



which in terms of the redefined density q becomes 



X J x-y J x-y L 

C23 C13 



Q{y) ^ 1 



COt23 



Cl3 



Acoti2 dy 
X — y 27rz 



27r/C23 + 4>2 



which coincides precisely with ( |3.4D as announced above! Thus the finite size corrections to 
the charge of any given configuration will indeed be equal to the field theoretical prediction, 
that is to the 1-loop shift around the classical solution. 



-17- 



5. Bosonic duality 

In this section we will explain some details behind the bosonic duality^ ( |2.12| ) mentioned 
in section |. There are two main steps to be considered. On the one hand we have to prove 
that for a set of K2 generic complex numbers U2 and Ki roots ui obeying the auxiliary 
Bethe equations ( p.3[ ) it is possible to write {t = <j)i — (j)2) 

2isin (t/2) Q2{u) = e"l'^Qx{u - i/2)Qi{u + i/2) - e-'^/'^Qi{u + i/2)Qi{u - i/2) , (5.1) 

and that, in doing so, we define the position of a new set of numbers ui. A priori this is not 
at all a trivial statement because we have a polynomial of degree K2 on the left whereas 
on the right hand side we have only K2 — Ki parameters to fix. However, as we will see, if 
Ki equations ( |2.3D are satisfied it is possible to write Q2{u) in this form. This will be the 
subject of the section 

Assuming ( |5.1| ) to be proved we can use this relation to show that in the original Bethe 
equations we can replace the roots ui by the new roots ui with the simultaneous exchange 
<Pi ^ <p2- Indeed if we evaluate the duality at n = U2j we find 

Qi{u2,j - i/2) ^ ^i(4>^-4>^) Qi{u2,j - i/2) 
Qi{u2,j + i/2) Qi{u2,j + i/2)' 



meaning that in the equation ( |2.4D for the U2 roots we can replace the roots ui by the dual 
roots ui provided we replace cpi ^ 4)2- Moreover if we take u = -uij it i/2 we will get 

^ Qi{ui+i) Q2{ui - i/2) 
Qi{ui-i) Q2{ui+i/2) ' 



which we recognize as equation (2^) with K2 — Ki roots ui in place of the Ki original 



roots ui and with 0i <-> (/>2. Finally evaluating {5A_) at u = uij it i/2 we will get the 
original equation ( p.3[ ) so that we see that it must be satisfied in order to equation ( |5.lD to 
be valid. 



In section |5.2| we will also see that the transfer matrices are invariant under the bosonic 

some 



duality accompanied by an appropriate reshuffling of the phases 4>a- In section 5.3 
curious examples of dual states will be given. 

5.1 Decomposition proof 

In this section we shall prove that one can always decompose Q2{u) as in (|5.1D and that 



this decomposition uniquely fixes the position of the new set of roots ui. In other words, 
let us show that we can set the polynomial 

P{u) = e+'iQi{u - i/2)Qi{u + i/2) - e'^^ Qi(u + i/2)Qi(n - i/2) - 2isin ^Q2{u) 

to zero through a unique choice of the dual roots Ui. 

^Bazhanov and Tsuboi also found some similar duality in the study of the deformed Uq{sl{l\l)). We 
thanks Z. Tsuboi for providing us the talk he gave at the "t9me rencontre entre physiciens theoriciens 
et mathmaticiens: Supersymmetry and Integrability" (http://www-irma.u-strasbg.fr/article383.html) and 
V.Kazakov who informed us of their work. It would be very interesting to connect both approaches. 
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Consider first the case Ki = 0. Then it is trivial to see that we can always find 
unique polynomial Qi = ""^^ + Yln=i ^nU^'^ such that 

e^'iQi{u + i/2)-e-'iQi{u-i/2) = 2i sin ^QsN • 

because this amounts to solving K2 linear equations for K2 coefficients a„ with non- 
degenerate triangular matrix. 

Next let us consider Ki < K2/2. First we choose Qi to satisfy Ki equations 



r 



Qi(u ) = 2ie 2 sm- \ ^ = Cp , p=l,...,Ki 

^ 2 Qi{u}p-i) 



these conditions will define Qi{u) up to a homogeneous solution proportional to 

Qi{u), 

where qi{u) is some polynomial of the degree K2 — 2Ki. Now from ( |2.3D we notice 
that with this choice of Qi we have 



P{ul + i/2) _ P{ul-i/2) 



, p=l,...,K3 



and thus 



where 



Q2{ul + i/2) Q2{ul-i/2) 

P{u) =Qi{u + i/2)Qi{u - i/2)p{u) 



■ T ■ T T 

p{u) = e^2q^(^u + i/2) — e~*2 — i/2) — 2i sin — q2{u) 



and q2 is a polynomial. Thus we are left to the same problem as above where Ki = 0. 
For completeness let us note that we can write q2 (u) explicitly in terms of the original 
roots ui and U2, 

, \ Q2{u) . 

Qi{u + i/2)Qi{u - i/2) 

where the last term is a simple collection of poles aX u = u^^i/2 whose residues are 
such that q2{u) is indeed a polynomial. 



We can see that the number of the solutions of ( p. 3D with Ki = K and Ki = K2 — K 
is the same (see for examples of states counting). Thus for each solution with 
Ki > K2/2 we can always find one dual solution with Ki < K2/2 and in this way 
we prove our statement for Ki > K2/2 

Finally let us stress the uniqueness of the Qi. If Ki > Ki we have nothing to show 
since we saw explicitly above how the bosonic duality constrains uniquely the dual 
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polynomial Qi. Let us then consider Ki < Ki and assume we have two different 
solutions Q\ and Qi- Then from the duality relation (5.1) for either solution we find 



e'iQiiu - i/2) (^Q\{u + i/2) - Ql{u + i/2) 
e-'^Qiiu + i/2) (qKu - i/2) - QHu - i/2) 

Evaluating this expression at tt = uij + i/2 we find that Q\{uij) — Qx(^ij) = so 
that Q\{ui) — Qi{ui) = Qi{u)h{u) and therefore 

e'ih{u + i/2) = e-'^h{u - i/2) 

which is clearly impossible for polynomial h{u) ~ for large u we can neglect the i/2's 
to obtain e*"^ = 1 thus leading to a contradiction. 

5.2 Transfer matrix invariance under the bosonic duality 

In this section we will examine the transformation properties of the transfer matrices under 
the bosonic duality. In Appendix A we consider this problem for the general su{N\M) 
group. For now let us just take Tn for su(l,2) from (^). Using (5T) we can express 
ratios of Qi's through Qi and Q2 so that 



Tu {u) 



+ e 



+ e" 



-t(p2 



+ 



2ism^e * 2^2(11 - | 
Q,{u+{)Q,{u + { 



Qi(u+|)Qi(n + |) 





2isin5e+^ig2(n + f) , +,,^.v^ , 4 



Qi{u + j 



We see that for t = (pi — (j)2 the terms with sin ^ cancel and we get the old expression for 
Tn with ui replaced by ui and (pi ^ (p2- 

This simple transformation property of the transfer matrices automatically implies 
that the Riemann surface defined by the algebraic equation ( |3.7D is untouched under the 
duality transformation (to all orders in L), so that the duality can cause at most some 
reshuffling of the sheets. However, as we will see in the next section, not necessarily the 
sheets as a whole are exchanged - this operation will be in general done in a piecewise 
manner. 

5.3 Examples 

In this section we will study some curious Bethe roots distributions for the twisted su(l, 2) 



^ij + 2 



spin chain described by the nested Bethe equations ( |2.3D and ( p.4[ ) and for the usual su{2) 
Heisenberg chain, 

L 

(5.2) 



Qi {uij + ■ 



Qi {uij 
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Using the first example we shall understand the importance of twists to stabilize big cuts 
of stacks like the ones depicted in figures 2a, 2b and explain how the stacks gets destroyed 
as we decrease the twists. 

We can dualize su{2) solutions of the twisted'' Heisenberg ring using the same duality 
(2.12) as before with Q2{u) — > u^- We will consider the dual solutions to the vacuum 



and to a 1-cut solution for the Heisenberg spin chain ( |5.2[ ) as a prototype of the curious 
solutions one would get. 

5.3.1 Big enough twists, small enough fiUings and zippers 

In the previous sections we saw that the introduction of twists in the NBA equations are 
needed to have a configuration with auxiliary roots ui close to some momentum carrying 
roots U2- In figure ^ we have two numerical solutions of the Bethe equations which are 
related by the bosonic duality. In either of them we see a configuration of Bethe roots 
with a simple cut with middle roots only (in blue) and a cut of stacks (containing blue and 
yellow roots). In this situation it is clearly reasonable to think of stacks as bound states 
of different types of roots and we see that they indeed condense into multicolor cuts. 

• * •••••••••••••••••••••••••••• jj_Q cccccoocooooooooo o o o 

^ 77^ 

Figure 6: The upper and the lower configuration of Bethe roots are dual to one another. Big blue 
dots are middle node roots U2, yellow dots are auxiliary roots ui. The formation of cuts of stacks 
is manifest for this situation where the twists are large (like 7r/2) and the filling fractions are small. 

We will examine what happens when we decrease the twists (or increase filling fractions, 
which is the same qualitatively). For simplicity we consider the configuration, dual to the 
simple one cut solution (K2 = K and Ki = 0) with no twist for the middle node roots, 
4'2 — (ps = 0) and some generic twist — (j)2 = t for the auxiliary roots. Bosonic duality 
will leave untouched middle node roots U2 and create K new axillary roots ni. 

In the upper left corner of figure we applied the duality for some big twist r = 4.6 
while in the bottom right corner of the same figure we have a configuration of Bethe roots 
with some small twist r = 0.2. In this latter case the auxiliary (yellow) roots clearly do 
not form stacks together with the middle node (blue) roots!, rather they form a bubble, 
containing the original cut of roots U2- 

To understand what happens in the scaling limit consider the position of 7123 = 1 



fluctuation, given by (4.3), which would be a small infinitesimal cut between p2 and p^. 
Clearly this probe cut would have no influence on the leading order algebraic curve for pi. 
In flgure |^ the position of this virtual fluctuation is marked by a red crossed dot. When the 
twist is big enough (and fllling fraction is small enough) the fluctuation is to the left from 
the cut. When we start decreasing the twist the fluctuation approaches the cut (upper 
right picture on fig 0) and at this point we have at the same time 

P2{Xn) -PsiXn) = 27r 

'''For zero twist the duality becomes degenerate and we will see below that it needs to be slightly modified. 
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T=4.7 T=1.7 



T=0.6 / \ T=0.2 



Figure 7: Disintegration of the stack configuration. Wlicn the twist is large (the top left corner) 
the auxiliary roots form bound states together with the middle node ones and constitute a cut of 
stacks. As we decrease the twist fluctuation 7123 = 1 (the red crossed dot) enters the cut of stacks 
(the top right corner) and subsequently partly disintegrate the cut of stacks forming some zipper 
like configuration (the bottom left corner) . At some very small value of the twist the configuration 
of Bethe roots bears no resemblance with a cut of stacks. 

and 

Pl{Xn) -PsiXn) = 27r , 

which implies pi — p2 = so that equation ( |3.13| ) becomes wrong at this point. When we 
continue decreasing the twist the fluctuation passes through the cut and becomes a ni2 = 
fluctuation. If we think of the fluctuation as being a small cut along the real axis we see 
that density becomes negative after crossing the cut: 

r, ^ flue MP2-P3) A(-pi-p2) flue 
< P23 = = = ~Pl2 

This means that two branch points of the infinitesimal cut should not be connected directly, 
but rather by some macroscopical curve with real positive density! This curves z{t) can be 
calculated from the equation p{z)dz £ M"*" or 

P^{z)-p,{z) ^ 
2m 

and the resulting curve is plotted in black on the two bottom pictures on the figure |^. 
This is very similar to what happens when a fluctuation passes through the 1 cut su{2) 
configuration [^4|. In the scaling limit the black curve corresponds to the cut connecting 
Pi and p2 like on the figure |^. 

At first sight these figures seem to be defying our previous results. Indeed we checked in 
the previous section that the transfer matrices themselves are invariant under the bosonic 
duality. Thus the algebraic curves obtained from (^]^) should be the same after and before 
duality and thus what one naturally expects is a simple interchange of Riemann sheets 
Pi ^ P2 under the duality transformation. What really happens is a bit more tricky. The 
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quasimomenta are indeed only exchanged but this exchange operation is done in a piecewise 
manner. That is, if we denote the new quasi-momenta by p"^"' and the old ones by pf''^ and 
if we denote the bubble in figure I by 7^ then we have 



new 

Pi 




outside TZ 
inside TZ 



new 
P2 




outside TZ 
inside TZ 



new 



^old 
Pi, 



where the border of the region TZ can be precisely determined in the scaling limit as 
explained above. 

5.3.2 Dualizing momentum carrying 
roots 

In this section we will consider an example 
of application of the bosonic duality to the 
Heisenberg mag net^. The duality ([2l^ ) 
can be applied to the roots ui obeying ( ^^ ) 
provided we replace Q2{u) . In fact 

if we want to consider strictly zero twist 
we need a new duality because that one 
is clearly degenerate in this limiting case. 
The proper modified expression is in this 
case 




Figure 8: In the scaling limit the algebraic 
curves for e*^^ are the same before the dual- 
ity (blue cut only) and after the duality (when 
the auxiliary roots are created). The duality 
causes interchange of the sheets outside the bub- 
ble, while keeping the order untouched inside. 
This follows from the need of a positive density 
for the "virtual" cut. In other words the duality 
is indeed only interchanging the sheets of the Rie- 
mann surface although it is interchanging them in 



i{ki-Ki)u^ = Qi{u-i/2)Qi{u + i/2)-Qi{u + i/2)Qi{u 



and the number of dual roots is now L — 
J^]*^^^^^fcon&ary to what happened with non-zero twists, here, the dual solution is not 
unique. Indeed if Ki > Ki we can as well use 



= aQi + Qi 



(5.3) 



All these solutions, parameterized by the constant a, have the same charges because 
the transfer matrix is invariant under this transformation - see appendix A. Notice that 
if initially we have a physical state with Ki < L/2 roots then all dual states ( |5.3| ) are 
unphysical with Ki > L/2 violating the half- filling condition. Still, it is interesting, at the 
level of Bethe equations, to understand how these solutions look like. First of all let us 
single out a particular Qi out of the various solutions to (^) so that 



1=0 

). We chose Qi 



(5.4) 



becomes well defined through (|5^ 

Consider for example the vacuum state for which Qi 
very large so that we can write 



Q\ to be the dual solution with Cg = 0. 

1. Let us first take a to be 



(5.5) 



*This section beneficed a lot from the insightful discussions with T. Bargheer and N. Beisert whom we 
should thank. 
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Figure 9: Three configurations of Betlie roots dual to tlie ferromagnetic vacuum of the untwisted 
Heiscnberg spin chain. For each physical solution (below half filling) of the Bethe equations there 
is a one parameter (a) family of dual unphysical solutions. To the left, a is large and the roots 
distribute themselves along a circle with radius Ra given by {RaL)^ — a. Decreasing a the circle 
will touch the fluctuations n — ±1. Similarly to the previous section the virtual infinitesimal cuts 
become macroscopical bubble cuts with cusps at the position of the fluctuations. Intersection points 
of the new cuts with the circle are connected by condensates, which are logarithmic cuts on the 
algebraic curve Q. 

We see for large a the dual roots will be on a circle of radius — . The corresponding 
configuration is present on the first picture on the figure ^. In this figure we also plotted a 
circle with this radius and one can see that the Bethe roots belong perfectly to the circle. 

Let us now understand this configuration from the algebraic curve point of view. The 
the quasi-momenta pi = —p2 = p = j^ — G, in the absence of Bethe roots, are simply given 
by p = Let us find the curves with positive densities and mode number n = 0. The 
density is given by p{x) = and we have to find the curves where p{x)dx is real. It is 
easy to see that the only possibility is the circle centered at the origin with an arbitrary 
radius. Prom the above arguments one can expect that for any a the roots will belong to 
some circle. However, we analysed only the curves with zero mode number and as we see on 
the figure ^ for smaller a's the circle develops four tails and two vertical lines. Along these 
vertical lines the roots are separated by i (for L ^ oo) forming the so called condensates 
or Bethe strings. The tails meet at the points where the virtual fiuctuation is and the 
corresponding curves are given by 

P^^d,z = ±1 
ni 

analogously to the previous section. In the last configuration on figure ^ the circle is 
completely absent. There are only two n = ±1 curves which, at the interceptions, become 
a 47r jump log condensate with the Bethe roots separated by i/2. 

We also built the dual configurations to the 1-cut solution (see figure |lO|) . The situation 
is similar to the vacuum, the only difference being that two tails (out of four) do not tend 
to touch each other, but rather end at the branch points of the initial cut. 
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(5.6) 
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Figure 10: Dual configuration to 1-cut solution. Similar to the previous example for the large a 
the dual roots are distributed along the big circle and cut (first picture). When the a decreases 
and the circle crosses the cut we have to choose another curve with the positive density (second 
and third pictures). 



6. The AdS/CFT Bethe equations and the semiclassical quantization of 
the superstring on AdS^ x 5"^ 

6.1 Introduction and notation 

The Beisert-Staudacher (BS) equations |2^] are a set of 7 asymptotic |Q Bethe equations 
(the rank of the symmetry group P5'C/(2, 2|4)) which are expected to describe the anoma- 
lous dimensions of = 4 SYM single trace operators with a large number of fields^ as well 
as the energy of the dual string states^'^. The perturbative gauge theory and the classical 
string regimes are interpolated by these equations through the t'Hooft coupling A. In 
based on an hypothesis for a natural extension for the quantum symmetry of the theory, 
Beisert found (up to a scalar factor) an S-matrix from which the BS equations would be 



derived. The scalar factor was then conjectured in |58, 59 1 from the string side - using 

and in from the gauge theory point of view - 

33i. From the gauge theory side these equations 
53, |6^, B^]. From the string theory point of 



the Janik's crossing relation |6C] 
based on several heuristic considerations 
were tested quite recently up to four loops 
view the scalar factor recently passed several nontrivial checks [67, 37, 68, ^ where several 
loops were probed at strong coupling. Also at strong coupling, the full structure of the BS 



equations was derived up to two loops in ||70|, |7l|] in a particular limit |7^] where the sigma 
model is drastically simplified. 

In this section we will check that the BS equations reproduce the l~loop shift around 
any classical string soliton solution with exponential precision in the large angular momen- 



^Tliese large traces can be tliougfi of as spin cliains and then tfie dilatation operator befiaves like a spin 
cliain Hamiltonian wfiidi turns out to be integrable In tliis way Betlie equations appear naturally 

from the gauge theory side. 

^°The existence of a finite gap description of the classical string motion ^ lead to the belief that 
these equations ought to be the continuous limit of some quantum string Bethe equations. In other words, 
the Riemann surfaces present therein should in fact be the condensation of a large number of Bethe roots. 
Inspired by these finite gap constructions these quantum equations were proposed shortly after |55l Eq| . 
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turn in the string state. To do so our computation is divided into two main steps. On the 
one hand we will compute the corrections to Bethe equations in the scaling limit. We 

will have to use the technology developed in the previous sections in order to understand 
precisely the several sources of corrections, the most subtle of all being the fine structure 
of the cuts of stacks which are generically present^^. At the end we will find out some 
integral equation corrected by a 1/\/A term. 

On the other hand we start from the algebraic curve description of the string classical 
motion The integral equations present in this finite gap formalism coincide with the 

scaling limit of the Bethe equations. Then we find how to correct this equations in such 
a way that they will now describe not only the classical motion but also the semi-classical 
quantization of the theory around any classical motion. For example we will find out how 
to modify the equations in such a way that they exhibit a very nontrivial property: the first 
finite corrections to any classical configurations equals the sum of quantum fluctuations 
around this same classical configuration. Then we show that, modified in this way, the 
integral equations coincide precisely with the scaling limit expansion of the BS equations 
with the HL phase |^9[ (up to some exponencially supressed wrapping effects, irrelevant 
for large angular momentum string states)! In this way we establish that, to this order in 
1 / \/A, the BS equations do provide the correct quantization of the system. 

These Bethe equations are a deformation of the equations ( p.lD through the introduc- 
tion of the map 

1 Attu _|_ 1 47r / i 

As explained in section |2| for superalgebras the choice of Bethe equations is not unique. In 
j25U four choices are presented. We need only to consider two of them^^, corresponding to 
the diagram in figure |l] or to the reflected path along the diagonal going from the lower 
left to the upper right corner. 

Moreover we consider a twisted version of these equations for the same reasons men- 
tioned in the previous sections. In |7C, 77] a similar kind of twists were introduced in the 
study of a set of deformations of = 4 SYM and of the dual sigma model. Our twists seem 
to be a simple change in boundary conditions via the introduction of a constant matrix 
like ( |2.2D . It would be interesting to see if they can also be given a deeper physical inter- 
pretation following the lines of these works. We should stress that the twists are used here 
as a technical tool which will simplify our analysis because, in particular, it allows us to 
deal with well defined stacks in a regime where the dualities are nothing but an exchange 
of Riemann sheets. We will explain in section 6.7 that we can then safely analytically 
continue the results to zero twist. 



^^In l?^, ^ the scaling limit of the 5(7(3) sector was considered. It would be interesting to use our 
treatment, including stacks, to compute explicitly the finite size corrections in this subsector following the 
lines of these papers. 

^^In psf we consider 771 =7^2 ~ Tj- 
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The BS equations then read 

iri4>i-iritl>2 



e 



^ir]<f)-i-iri<i>4. 



e 



7=1 ^3,fc - ^^2,i - § a;3,fe - xlj ' 



1 1 / + II II II ~i 1 / + ' 



^ir]<f)Q-ir}<f)7 



-j-j- U^,k — Uqj 2 ~ ^4,j 

j=l ^5, A; ~ "^GJ ~ 2 j=i ■^5,k ~ X^j 

TT ^6,fc - Uii,j - ^ TT ^e.fc - ^5,j + 2 TT ^6,fc - '"Tj + 2 



7=1 ^^7,fc - U6J - f ^^J^ 1 - l/xj^kX^j 



In fact, in order for the fermionic duality |^] (which we will review below) to exist, the 
twists must not be completely independent but rather 



^ 1 x^ 

01 - <?!'2 + V - log ^ = (t>3-(P4 , 
(/)7 - 08 + - log ^ = 05 - 06 ■ (6.2) 



The energy (the anomalous dimension) can then be read from 



j=l \^4,j 



To describe classical solutions (and to semi-classically quantize them) we should consider 
the scaling limit where 

~ U ~ Ka ~ L ^ 1 . 

In this limit we have 

x±'-a{x) + o(j 
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where 

, , 47r 

It is then useful to introduce the resolvents^^ 

^ U- Ua,j 



P5 = 2 \ + ^ (-^5 + ^4 - ^7) + <^5 



(6.4) 



and build with them eight quasi-momenta {^J = L/\/A) 

P. = " '-""^^^^ + V i-H, + H, + H,- H,) + 

P4 = ''"-^^^^ + - + + 04 
x^-l 

P7 = - ^"-^^ " ^-"^.^ ''""^ , (-Fe + + - Fe) + 0. 

where G4(x) = — Yl^=o Qn+ix^. We can also write 

Then, to leading order, these quasi-momenta define an eight-sheet Riemann surface and 
the BS equations read simply j/j — = 27rnij in each of the cuts Cij uniting pi and pj. 
Finally, in this section we will use 

p'i ~ P'' Pi ~ P' 
cotij = a{x)^—^ — - cot -^—^ — - 

which is similar (but should not be confused) with (|3 



note that 

Faix) = Gaix) + Gaix) = Ha{x) + Ha{x) . 
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Figure 11: The several physical fluctuations in the string Bethe ansatz. The 16 elementary physical 
excitations are the stacks (bound states) containing the middle node root. From the left to the 
right we have four 5'^ fluctuations, four AdS5 modes and eight fermionic excitations. The bosonic 
(fermionic) stacks contain an even (odd) number of fermionic roots represented by a cross in the 
psu(2,2|4) Dynkin diagram in the left. 



6.2 Middle node anomaly 

In this section we will expand BS equations in the scaling limit for the roots belonging to 
a cut containing middle node roots X4 only. We do not assume that all the others cuts are 
of the same type, rather they can be cuts of stacks of several sizes. In the section we 
will generalize the results obtained in this section to an arbitrary cut, assuming, as in the 
previous section, that the cuts are small enough and twists are not zero so that stacks are 



stable. We will discus in section 6.7 what happens when one takes all twists to zero. 

To leading order, the middle node equation ( |6.1D can be simply written as ^4—^5 = 27rn 
while at 1-loop the first product in the r.h.s. of (|6.1[) corrects this equation due to 



K4 , 



2 ^4{x) + a{x)TTp'{x) COt{TTp{x)) (6.5) 



where p{x) = Expansion of the remaining terms in (|6.1| ) will not lead to the appearance 
of such anomaly like terms since the roots of another types are separated by ~ 1 from x^^k- 
Thus we have simply 

27rn = ^4 — ^5 — r/ a(x)7r/)'(x) cot(7r/)(x)) , a; G C45 (6.6) 

In the next sections we will use dualities of the BS equations to get some extra information 
about cuts of stacks and generalize the above equation to any possible type of cut. To 
achieve this we shall recast this equation in terms of the middle node roots X4 only. 

6.3 Dualities in the string Bethe ansatz 

Obviously, the behavior of the Bethe roots will be as described in section |^ for a simpler 
example of a sti(l, 2) spin chain, that is, we will have simple cuts made out of X4 roots only 
and also cuts of stacks with X2,X3 and X4 roots for example. Consider such cut of stacks. 



Clearly, to be able to write the middle node equation (|6J|) or (|6.6| ) we need to compute 
the density mismatches p2 — P3 and ps — p^ which are 1-loop contributions we must take 
into account if we want to write an integral equation for the middle node equation in terms 
of the density /34 of momentum carrying roots only. In this section we shall analyze the 
dualities present in the BS Bethe equations. By analyzing them in the scaling limit we will 
then be able to derive the desired density mismatches. 

6.3.1 Fermionic duality in scaling limit 



In [25 1 it was shown that the BS equations obey a very important fermionic duality. Since 
we chose to work with a subset of the possible Bethe equations, that is the ones with 
f]! = f]2 = V present in |2^, we should apply the duality present below not only to the 
fermionic roots xi and X3 (as described below) but also to the Bethe roots X5 and X7. 
Obviously the duality for X5 and xj is exactly the same as for xi and X3 and so we will 
focus simply on the latter while keeping implicit that we always dualize all the fermionic 
roots at the same time. 

We construct the polynomial {t = rj (04 — cps)) 

Ka K2 
P{x) = e+*^ \{{x - x^.) JJ(x - xl^){x - 1/X2j) 

Ka K2 

- e"*^ - ^r,.-) - - (6.7) 

of degree + 2K2 which clearly admits x = x^j andx = 1/xij as K3 + Ki zeros^^ . The 
remaining K/j^ + 2K2 — — Ki roots are denoted by x^j or 1/xij depending on whether 
they are outside or inside the unit circle respectively, 

ifi Ki K3 K3 

P(x) = 2isin(T/2) J](x - 1/xi,,) Y[{x - l/xi,,) JJix - xgj) - ^3,i) (6.8) 
j=i j=i j=i j=i 

Then we can replace the roots xij,X3j by the roots xij,X3j in the BS equations provided 
we change the grading r] ^ —rj and interchange the twists (pi <-> (/>2 and (/>3 <-> (/>4. In fact, 
since we should also dualize the remaining fermionic roots, we should also change (j)^ <-> (pQ 
and (pi (ps and replace the remaining fermionic roots X5 and X7. 

Since to the leading order x^ ~ x each root will belong to a stack which must always 
contain a momentum carrying root X4. We have therefore Ki = K2 — Ki and K^^ = 
K2 + K/i — K3. Thus we label the Bethe roots as 

xi,j = X4j - eij , j = 1, . . . ,Ki 

xi,j = - h,j , j = ■ ■ ■ ,Ki 

X2,j = X4J - e2,j , j = 1, . . . ,K2 

X3,j = Xij - eaj , j = 1, . . . , K3 

X3,j = Xi^j+Kz - h,j , j = ■ ■ ■ ,K3 



14, 



we also have 1/xi has zeros because, due to (6.2), the equation for xij is the same as the equation for 



X3J if we replace x^j by 1/xij. This is why the restriction (3.2) of the twists is so important 
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with e ~ 1/VA. Dividing Q and Q by Uf=ii^ - ^'ij) Uflii^ - X4,j){x - l/x^j) we 
have 



a-2 j J- / 2,j -iL "TT ■^4,j "TT ■^2,j ' 2,j 



= 2i sin(T/2) TT ^^iZ^ TT rr TT (6.9) 

J=l ^ ~ V2;4,j fj-^ X - l/xA,K,+j fj^ X - X4J j-J^ X - XA^K^+j 

In this form it is easy to expand the duahty relation in powers of \/\/~\. By expanding all 
factors in (|6.9D such as 



TT ^ = exp V log ^ ~ exp T7;G2{x) + V — — 

fJ-^X-XiJ yj^^ X-Xijj y 2 ^X-X2j 

we find 



2 / \ ^-^ X — X3 ^ a; — X3 X — X2 

V ^x-l/xi ^x-l/x2 



Then, similarly to what we had in section 3.2 for the bosonic duality, we notice that 
a{x)dJ^^ + Y.^-T.^)=H^ + H3-H^-H^^ 

Y'^^ X — X3 X — X3 X — X2 J 

with a similar expression for the argument of the second exponential. Thus finally we get 

{Hi + H2-H;- H~^) + {H2 -Hi- H~^ = -cot34 , 

or alternatively, using the x ^ 1/x symmetry transformation properties of the quasi- 
momenta, 

(^4 + H2-H3- ^3) + {H2 -Hi-Hi) = - coti2 . 

From this expressions we can deduce several properties of the density mismatches we wanted 
to obtain. For example, if we compute the discontinuity of ( |6.3.1 ) at a cut containing roots 
xi, that is in a large cut of stacks Ci^i>4, we immediately get 

Acoti2 , - 

Pi- P2 = , X G Ci,i>4 . (6.10) 

Proceeding in a similar way we find 

A COt34 ^ ^ . „ , ^ . 

P^- Pa = ^— , X G C3,i>4 , (6.11) 

P3 - Pa = P2 - P3 , X e Ci,i>4 U C2,i>4 • (6-12) 

Let us now show that in the scaling limit the fermionic duality corresponds just to the 
exchange of the sheets {pi} of the Riemann surface. For illustration let us pick pi and see 
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Figure 12: Action of tlie duality on a long stack. By successively applying the fermionic and the 
bosonic dualities duality we can reduce the size of any large cut. One should not forget to change 
the sign of the grading 77 after applying the fermionic duality. 



how it transforms under the duality. By definition the fermionic duality corresponds to the 
replacement tj —r], Hi — > Hi^H^ and (/>2, 03 </'4, so that 



Pi 



2ttJx — 6r]-lQl + 6ri,+lQ2X 



In the same way we get 



T] {-Hi - Hr^ + Hi) +(l)2=P2+r] COtl2 



P2^Pl+'n COtl2 , Pa^PA-V COt34 , P4 ^ P3 - V COtu , 

and since cot^ ~ 1/\/A we see that to the leading order the duality indeed just exchanges 
the sheets. 

6.3.2 Bosonic duality in scaling limit 

The bosonic nodes of the BS equations are precisely as in the usual Bethe ansatz discussed in 
the first sections so that we can just briefly mention the results. The duality (r = r]{(p2—4'3)) 



e+'2Q^{u - i/2)Q2{u + i/2) - e~'2Q^{u + i/2)Q2{u - i/2) = 2ism -Qi{u)Q3{u) 



leads to 



which implies 



{Hi + HS-H2- H-^) + {Hi + H3-H2- H^) = cot23 



, A COt23 ^ ^ 

P2- P3 = H TT— , X £ C2,j>4 

ZTTl 



(6.13) 



As we already discussed in section ^ the bosonic duality also amounts to an exchange 
of Riemann sheets. Indeed, under the replacement H2 H2 and 4>2 ^ (ps^ we find 

P2^P3-r] COt23 , P3^ P2 + r] COt23 

which again, to the leading order in \/A, is just the exchange of the sheets of the curve. 
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C2,i 


^3,1 


27ri(/9i 


- P2) 


— A COtl2 






2'Ki{p2 


-Ps) 


—A cot 13 


+A COt23 




2vri(/93 


-P4) 


+A cot 14 


— A COt24 


— A cot 34 



Table 1: Densities missmatches 



6.3.3 Dualities and the missing mismatches 

Using bosonic and fermionic dualities separately we already got some information about 
the several possible mismatches of the densities inside the stack. To compute the missing 
mismatches we have to use both dualities together. For example suppose we want to 
compute p3 — /94 in a cut Ci^i>4. We start by one such large cut of stacks (see figure 12a) 
and we apply the fermionic duality to this configuration so that we obtain a smaller cut as 
depicted in figure 12b. For this configuration we can use ( |6.3.2| ) to get 

A cot 14 



However, from ( 6.12 ), this is also equal to the mismatch we wanted to compute, that is 

A cot 14 

/03 - P4 = H — TT— ' ^ £ Ci,j>4 . 
Zttz 

To compute the last mismatch we apply the bosonic duality to get a yet smaller cut as in 
figure |l^ c for which we use ( |6.11 ) to get 

A cot 13 

Again, from ( |6.12 ), we can revert this result into a mismatch for the configuration before 
duality, that is 

A cot 13 

P2- P3 = TT— 1 ^ £ L-l,j>4 . 

Zttz 

Let us then summarize all densities mismatches in table [l|. 
6.4 Integral equation 

In this section we shall recast equation (^) or 

V ; \-2^4-Hs-H5-Hi-H7 = 27rn+7/(/)4-?/(/)5-cot45 (6.14) 

— 1 

in terms of the density p4{x) of the middle roots X4. To do so we only need to replace the 
several densities by the middle node density P4{x) using the several density mismatches 
presented in table |l[ Defining 
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we can then rewrite equation ( 6.14| ) in terms of the middle node roots only, 



V ^^'^9 ^' — ^ 2 -^45 + -f^is + H^s — 2His — — H4S 

— 1 

= 27m + r/<^4 - V<P5 - cot45 + Yl ^ +^3)+ E ^^li + ^^') (6.15) 



l<i<4 l<i<4 
5<i<8 5<j<8 



where x S C45 and 



The several dualities amount to an exchange of Riemann sheets so that the cuts Cij Ci'j' 
with the subscripts in Hij changing accordingly. The middle roots X4 are never touched in 
the process. Moreover to leading order pi <-> p^/ and thus the r.h.s. of ( |6.15| ) is also trivially 
changed under the dualities. Therefore, as in section |3| (see (3^) and (|3.5|)), we can now 



trivially write the corrected equation when x belongs to any possible type of cut of stacks 
by applying the several dualities to equation (6.15). 



6.5 Fluctuations 



In this section we shall find the integral equation ( 6.15| ) from the field theoretical point of 



view like we did in section and in appendix B. That is, we will find what the corrections 
to the classical (leading order) equations 

V :; ^ 2^4 - H3- H5- Hi - Hj = 27rn + rjcp^ - r]<p5 , (6.16) 

— 1 

should be in order to describe properly the semi-classical quantization of the string (and not 
only the classical limit). We will find that this construction leads precisely to the integral 
equation ( |6.15| ) thus showing that the BS nested Bethe ansatz equations do reproduce the 
1-loop shift around any (stable) classical solution with exponential precision (in some large 
charge of the classical solution). This section is very similar to section ^ and to Appendix 
B and thus we will often omit lengthy but straightforward intermediate steps. We assume 
i = 1, . . . , 4 and j = 5, . . . , 8 in all sums. 

As in ( [4.6| ) and ( [7.6D , we add ^(—1)^ of a virtual excitation for each possible mode 
number n and polarization ij to each quasi-momenta. Notice that for this super-symmetric 
model the fiuctuations can also be fermionic and indeed the grading (—1)^ equals +1 (—1) 



for bosonic (fermionic) fluctuations, see figure 11, as usual for bosonic (fermionic) harmonic 
oscillators. 

We denote p = po + ^p where po is the leading density, solution of the leading (classical) 



equation ( 6.16 ), while p obeys the corrected (semi-classical) equation. For example, if we 



consider x G ^4,5, the starting point should be (see for a similar analysis) 

-2x6n-i6Qi ^ ^ f a{x) 6p{y) ^ f a{x) 5p{y) ^ f a{x) 5p{y) 

x^ -1 J a{y) x-y J a{y) x-y J a{y) x-y 

C45 C15 C48 
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a{l/x) 5p{y) f a{l/x) 6p{y) f a{l/x) 6p{y) 



a{y) l/x-y J a{y) l/x-y J a{y) l/x-y 

Clg Ci5 C48 



^ 1 



a(x) 



a{l/x) 



n=-N 

Then, by construction, the charges 

piy) 

ay ^ ^ 



j<4 j>5 



a(l/x 

j>5 ^/""^ ~ 



(6.17) 



Qr 



c y 



'-^dy+yy{-if^^^ 



p{y) 



c y 



dy + Y. 



cotjj dy 
y^ 2TTi 



(6.18) 



win take the l/^/X corrected values. It is clear that, as before, we do not include the 
new virtual excitations in the density p{x). Similarly to ( [4.12| ) and ( [7.1l| ), if we want the 
charges to have the standard form 



Qr 



Qiy) 



dy 



we must redefine the density as 



^( E (-l)'''''Acoti,,+ y (-l)^..'Acot,y j 



Q = P + 



Now we want to go back to the integral equation ( |6.17|) and rewrite it using the density 
6q = Q — pq. For example, for x £ C45, 

2 I a(a^) Sp{y) ^ f a{x) 6p{y) ^ f a{x) 5p[y) 



a{y)x-y J a{y)x-y J a{y) x - y 
C45 C15 C48 



N 



^^^^^ ^2^'^^ ^^^^^ 



n=-N 

2 y q(x) bQ[y) ^ r a{x) 6Q{y) ^ f a{x) 60{y) 



4:j 



a{y)x-y J a{y)x-y J a{y) x - y 

C45 Cl5 C48 

+ COt45-E(^+^)-^E(^+^/+^+^ 



where the identity 



(-1)^^ cot4, = - E +^) - E C^^^ +^ ' 



where i = i, 1 = 4, 2 = 3, is being used. Now, when x G Cis, we will get 
2 I Q(^) ^Pjy) _^ f "(a;) Sp{y) , f oi{x) 6p{y) 



a{y) x-y J a{y) x-y J a{y) x-y 
C18 C15 C48 
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N 



+ E 

n=-N 

a{x) Sgiy) 



+ 



a(x) 6Q{y) f a{x) 6g{y) 



+ 



a{y)x-y J a{y)x-y J a{y) x - y 

Cl8 Cl5 C48 



^ /^ri* -u T-J's t^J 

9 ~^8i -^H 



Finally we can use the x to 1/x symmetry to translate last equality into one for x £ C45. 
Subtracting it from the previous equation we see that the 1/\/A corrected equation will 
correspond to adding 

-cot45 + j;(2f +2^+2Ji+2^) 

to the r.h.s. of ( |6.16| ) thus obtaining, after the identification g = ^34, precisely the finite 
size corrected equation ( |6.15| ) obtained from the NBA point of view! 

6.6 The unit circle and the Hernandez-Lopez phase 

In the last section we showed that the one loop shift as a sum of all fluctuation energies 
(or others local charges) perfectly matches the finite size corrections in the NBA equations. 
However we systematically dropped the contours around the unit circle. 



For example, when we blow the contour in the last term of (6.18), we also get some 
contribution from the singularities inside the unit circle. That is we will have an extra 



contribution to the charges given by an integral over the unit circle. Also, take (6.17) for 
instance. To pass to the r.h.s we transformed the collections of poles into integrals over 
the excitation points and then we blew the contour which became a collection of contours 
on the several existing cuts. Again we dropped the contribution from the integrals over 
the unit circle which would lead to an extra term in the r.h.s. of ( |6.15 ). In our 



previous paper |32] we showed^^ that this extra contribution matches precisely the extra 
contribution coming from the Hernandez-Lopez phase in the NBA! 



However, as we explained in |32|, in order to obtain precisely the HL phase a precise 
prescription for the labeling of the mode numbers of the fluctuations must be given. 

Moreover, in we assumed that everywhere we can replace cot (^EiMzEiM^ by 

i sign(Im x) with exponential precision in This is reasonable for generic points in the 
unit circle, where the imaginary part of pi{x) — Pj{x) is large, but one has to carefully 
analyze the neighbourhood of the real axis, where this imaginary part vanishes. 
Let us consider these two subtle points in greater detail. 

6.6.1 A mode number prescription 



As we emphasized in [31| if we number the fluctuation charges differently we might 
obtain different results for the 1-loop shift, that is for the graded sums of these fluctuation 



Recently the HL phase was also found jT^] in the study of the open string scattering of giant magnons 
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charges. Thus a precise prescription for the labehng of the quantum fluctuations is crucial. 



In the appendix A of 1 32 ] we found out that the contribution of the integrals of the previous 
section does reproduce the HL phase provided we number the quantum fluctuations located 
at Xn according to 

Pii^li) - Pj{xi) =2TT{n-mi + rrij) 

with some specific choice of rm . Moreover we also showed that for the same choice of rm the 
contribution to the charges coming from the above mentioned integrals over the unit circle 
is zero. Using the x to 1/x symmetries following from the definition of the quasi-momenta 



(6.4) plus the restriction ( |6.2[ ) on the twists, we can redo the computation in the Appendix 
A of to find that the condition on the rrii now reads 

(m2 + — mi — 711,4) (w-5 + ^8 ~ ^6 ~ f^i) = 



so that, in particular, mj = does the job nicely. We see that, with the introduction of these 
twists and subsequent redefinition of the quasimomenta, the prescription for the labeling of 
the excitations becomes absolutely natural and algebraic curve friendly [^]. This answers 



the question raised in |32| concerning the naturalness of the presciption needed to obtain 
the HL phase 



see appendix A in |32]. 
6.6.2 Unit circle contribution 

Let us now us focus on the vicinity of x = 1 where we have the following expansion of the 
quasi-momenta 

2 ~ x-l 

where Pij is usually of order L/^/X (and should be so for the asymptotical BAE to be 

where is some large 



valid). We will consider the circle with radius x^^_^_i/2 



1 + 



cutoff in the sum of fluctuations (3.17). We want to estimate 



TTNl3i. 



a{x)f{x) 



cot 



Pi - Pj 



+ i sign(Im x) 



{p'i -p'j)dx. 



This integral is dominated for x ~ ±1 and can be performed by saddle point. The contri- 
bution for X ~ 1 is 



/ 



a{x)f{x) 



cot 



Pi - Pj 



+ i sign(Im x) 



{p'i-p'Mx 



i^'/(l) 
6/3i,VA 



+ 



which is zero under the sum over all polarizations. For example 

(_ 1)^45 (-1)^35 



'45 



'35 



Thus we can indeed drop the cot's when integrating over the unit circle and thus we finally 
conclude that the one loop shift to any local charge computed from the BS equations with 
the Hernandez-Lopez phase is indeed given by the sum of fluctuations as predicted by field 
theoretical arguments. 
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6.7 Zero twist and large fillings via analytical continuation 



Although we always assumed the twists to be sufficiently large and the fillings to be suf- 
ficiently small we can always analytically continue the results towards zero twists or large 
filling fractions. Let us briefly explain why. In the scaling limit, for large twists, the bosonic 
duality we introduced amounts to a simple exchange of sheets in some Riemann surface, 
Pa{x) Pb{x)- As we saw in section 5.3 what happens when the twists start to become 
very small is that the quasi-momenta are still simply exchanged but in a piecewise manner, 
that is, we can always split the complex planes in some finite number of regions where the 
bosonic duality simply means Pa{x) Pb{x)- Thus, from the e*^ algebraic curve point of 
view nothing special occurs for what analyticity is concerned and therefore we can safely 
analytically continue our findings to any value of the twists. Exactly the same analysis 
holds for the filling fractions. Moreover, for the usual Bethe system, we defined a set of 
quasi-momenta, which constitute an algebraic curve to any order in 1/L, and therefore we 
don't expect analyticity to break down at any order in 1/L. 

We also preformed a high precision numerical check concluding that there is no sin- 
gularity when the configuration of the Bethe roots is affected by this partial reshuffling of 
the sheets and that finite size corrections are still related to the same sum of fluctuations, 
which are analytical functions w.r.t. the twists. 



7. Conclusions 

In this paper we studied generic nested Bethe ansatz (NBA) equations, the corresponding 
scaling limit and its leading finite size corrections. Let us summarize briefiy our main 
results 

• We found out that the introduction of some extra phases, called twists, are crucial 
for the formation of bound states of roots of different types, called in the literature 
by stacks [|^]. Strictly speaking without these twists the stacks do not exist. See 
sections I and 

• We understood how to use the bosonic duality between various systems of Bethe 
roots which is present even in the absence of any fermionic symmetry. In the scaling 
limit we showed that this duality amounts to a reshuffling of Riemann sheets of the 
algebraic curve formed by the condensation of Bethe roots. See sections § and ^. 

• We explained how to write down the integral equation describing the leading finite 
size corrections around generic NBA's for (super) spin chains by using the transfer 



matrices for (super) group along with some TQ relations. See section 3.6 



We provided an alternative derivation of this integral equation using an independent 
path, namely using the dualities present in the Bethe equations allowing one to get 
rid of the several stacks and reduce the size of any cut by successive application of 



several dualities. See section 3.2. 
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We obtained the integral equation describing the finite size corrections to the Beisert- 
Staudacher equations |25| with the Hernandez-Lopez phase |11, 39 1 in the scahng hmit 
(to do so we were forced to use the duahty approach because at present the psu{2, 2|4) 
transfer matrices for this Bethe ansatz are not known^^). See section ^. 



In the scahng hmit Beisert-Staudacher equations [25| describe the classical motion of 
the superstring on AdS^ x through the finite gap curves of 1^]. Thus the integral 
equation we found should reproduce the 1-loop shift for all the charges around any 
classical string motion and this is obviously a very nontrivial check of the validity of 
the BS equations. We show that this equation indeed mimics the presence of a sea 
of virtual particles thus proving this general statement. See section I 
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Appendix A: Transfer matrix invariance and the bosonic duality for SU{K\M) 
supergroups 



In this section we review the formalism of [24| which allows one to derive the transfer ma- 



trices of usual (super) spin chains in any representation. We will use this general formalism 
to prove the invariance under the bosonic dualities of all possible transfer matrices one can 



build. The transfer matrices presented in section S.C can be obtained trivially using this 
formalism^^. 

As mentioned in section ^, for the standard SU{K\M) super spin chains (based on the 
standard i?-matrix R(u) = u + iV with V the super permutation) we can find the (twisted) 
transfer matrix eigenvalues for the single column young tableau with a boxes through the 



non- commutative generating functions [p4| , 4C] 



^^See section 6 in for some attempts to fill this gap. 

^'^We should mention that the transfer matrices in section 5.6 are not exactly the same we have in this 



Appendix but can be obtained from these via a trivial rescaling in u which obviously does not spoil the 
invariance of these objects. 
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where 7 is a path starting from (M, K) and finishing at (0, 0) (always approaching this 
point with each step) in a rectangular lattice of size M x i^T as in figure ||^^, x = (m, k) is 
point in this path and n = (0, —1) or (—1, 0) is the unit vector looking along the next step 
of the path. Each path describes in this way a possible Dynkin diagram of the SU{K\M) 
super group with corners denoting fermionic nodes and straight lines bosonic ones, see 
figure g. Finally, 

y^i ^ gi</,, Qk,m{u + i{m-k- l)/2) Qk-i,miu + i{m - k + 2)/2) _ ^-g^ 

M,io,~i^ ' Q^^^(u + i{m-k + l)/2) Qk-i,m{u + i{m-k + 0)/2) 

1 _ fj^^ Qk,m-i{u + i{m-k- 2)/2) Qk,m{u + i{m - k + l)/2) _ -^^ 



^ V Qk,m-iiu + iim-k + 0)/2)Qk,Uu + i{m-k-l)/2' 

where Qk,m is the Baxter polynomial for the roots of the corresponding node^^ and {(pk, ^m} 
are twists introduced in the transfer matrix pO| . Let us then consider a bosonic node like 
the one in the middle of figure || (the vertical bosonic node is treated in the same fashion). 
If the position of this node on the M x K lattice is given by (m, k) then it is obvious 
that the only combination containing Qm,k in the right hand side of ( [7.1[ ) comes from the 
product of (_, o)(^i)t>7„^+i^^) (_, o)(n) which reads 

Qk,m+i{u + i{m-k + 2)/2) Qk,m-i{u + i{m - k - 2)/2) ^ ^^■q^_ 
Qk,m+i{u + i{m - k + 0)/2) Qk,m-i{u + i{m -k + 0)/2) 
ii^^+i Qk,m{u + i{m-k- l)/2) Qk,m+i{u + i{m-k + 2)/2) 

(7.2) 



Qk,m{u + i{m - k + l)/2) Qk,m+i{u + i{in - k + 0)/2) 
I gi^^ Qk,m-i{u + i{m-k + 0)/2) Qfc,m(^ + i{m - k + 3)/2) \ ^-g^ 
^ Qk,m-i{u + i{m-k + 2)/2) gfc,„(u + i(m - + l)/2) J ^ 



So, if we want to study the bosonic duality on the node {k, m) and its relation with 
the invariance of several transfer matrices we need to study the last two lines of this 
expression. For simplicity let us shift omit the subscript k in the Baxter polynomials 
Qk,m-i,Qk,miQk,m+i and define the reduced transfer matrix as 

Qmiu- i) Qm+liu + if^) , J^^ Qm-l{u - i/2) Qmju + j) 



Qmiu) Qm+l{u-i/2) Qm-l{u + i/2) Qm{u) 

(7.3) 

Notice that the absence of poles at the zeros of Qm yields precisely the Bethe equations 
for this auxiliary node. 

Bosonic duality ^ Transfer matrices invariance 

Thus, to check the invariance of the transfer matrices in all representations it suffices to 
verify that the reduced transfer matrix t{u,ipm,V'm+i) is invariant under Lfm <-> fm+i and 

Notice that the path goes in opposite direction compared to the labelhng a of the Baxter polynomial 
Qa used before. In the notation of this section Qk,m corresponds to the node is at position (m,k) in this 
lattice. 

^^Qo,o is normalized to 1. If we are considering a spin in the representation where the first Dynkin node 
has a nonzero Dynkin label then Qm,k will play the role of the potential term. In general the situation is 
more complicated, see . In any case we are mainly interested in the dualization of roots which are not 
momentum carrying thus we need not care about such matters. 
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Qm Qm where 



2,sin i^flIl±l_£lIL^j Q^_,{u)Qm+i{u) = (7.4) 
e' 3 - i/2)Qm{u + i/2) - e"* 2 Qm(M + i/2)Q^{u - i/2) . 



which can be easily verified. If suffices to replace, in t{u, (pm, fm+i) in (|7.3|) , 

Qm{u-i) ^ ^-i(^^^^-^^) Qm(.U-i) 
Qmiu) Qm{u) 

frn+l-frn . f - Qm-l{u + 1 / 2)Q m+l{u + t /2) 

+2ie 2 sm ' ' 



2 / Qm{u)Qm{u) 

Qm{u + i) ^ g+i(^„+i-^„)OmOi + i) 
QmH Qm{u) 

'Pm+l-frn . - V'm\ Qm-ll"" - «/2)Qm+l(^^ - V2) 

— z«e 2 sm ' ' 



2 / Qm{u)Qm{u) 

which are obvious consequences of the bosonic duality. 
Transfer matrix invariance Bosonic duality 

On the other hand suppose we have two solutions of Bethe equations, one of them character- 
ized by the Baxter polynomials {. . . , Qm-i, Qm-, Qm+i-, • • •} with twists {. . . , ^m-, V'm+i, ■ ■ ■ 
and another with {. . . , Qm-i,Qm, Qm+i, • • •} with twists {. . . , ^fm+i, ^m, . . .} for which the 
transfer matrices are the same, that is 

t{u,(fm,fm+l)=i{u,(pm+l,fm)- (7.5) 



Then we can show that these two solutions are related by the bosonic duality (7.4). Indeed 
if we build the Wronskian^*^ like object 

_ . >Pm+i-'Pm Q^{u - i/2)Qm{u + i/2) _. ^m+i-^m Q ^^{u + i/2)Qm{u - i/2) 

yy ( Hi ) — Q, 2 — ^ 2 ^ 

Qm—l 

we can easily check that 

W{u + i/2) - W{u-i/2) = 

Qm-i{u-i/2)Qm+i{u + i/2) 

Since by definition W{u) is a rational function this means it must be a constant. Thus 
if ifm 7^ H^m+i we must have Km + Km = Km + Km+i and the value of W can be read 
from the large u behavior. In this way we obtain precisely the bosonic duality (|7.4| ). If 
'■Pm = Pm+i then we see that Km + Km = Km + Km+i + 1 and we will obtain a different 
value for the constant W which will correspond to the untwisted bosonic duality described 



in section 5.3.2 



2°We would like to thank A.Zabrodin and V.Kazakov for sugesting this nice interpertation for the bosonic 
duality 
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Appendix B: Fluctuations for su{n) spin chains 



In this Appendix we consider a su{n) NBA with the Dynkin labels Va being +1 for a 
particular a only (the generalization is obvious). This example is obviously more general 
than that considered in section O and can be a useful warmup for section where we 
find the integral equation describing the AdS^ x 1-loop quantization. For the spin 
chain su{n) NBA, in the classical limit, we will have n quasi-momenta each one above or 
below each of the n—1 Dynkin nodes^"*^. We label these quasi-momenta by pi (pj) with i, i' 
taking positive (negative) values for quasi-momenta above (below) the node for which 
Va ^ 0. Then let us mention how the equations in the previous section are generalized. 
We consider a middle node cut Ci-i. The analogue of equation ([4. 61) is now 



x-y 



+ 



x-y 



+ 



N 

E 



-N 



1 

2L 



■Vt + E 



and the charges (|3), (gj), (|4|) and ( pi)) b ecome 



X — X 



(7.6) 



Qr 



p{y) 



dy 



1 1 _ 1 \^ 1 



+ 



-E 

2L ^ 



cot 



y 



13 dy_ 
2m 2L 



-E 



cotjj dy 
i y^' 2-111 

cot ij dy 



y 



27ri 



-y 

2L ^ 



cotjj' dy 



1 



y"- 2m 2L 



-y 



cotjj/ dy 



2L Jc V 
so that the natural definition of the dressed density becomes now 



dy 



2m 



2m 



(7.7) 
(7.8) 
(7.9) 
(7.10) 

(7.11) 



Next step is to rewrite the integral equation ( [7.6D in terms of this new density. We proceed 
exactly as in ( 4.13| ), ( [4.14| ) using now 

cotfci(y) dy 

' Ci'i 




coti 



^ Ci 



y 27ri ' 



which is the analog of (|3.13| ) for this su{n) setup, so that at the end we obtain the following 
equation 



E 



Ci 



My) _^y- 
X — y ^ 



h- COtl_l 

X — y L 



E 



A coti^i -|-A cotj^_i dy 
X — y 2TTi 



(7.12) 



^'^See figure for an example of sucfi pattern for a super group wliicli clearly resembles su{S). 

^^as in the previous section, we are ignoring the regularization of the charges coming from the contribution 
of the contour around the origin which would appear in the second line from opening the contours around 
the excitation points a;^. 
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for Sq = g — go where go obeys the leading order equation 



M+y f ^ = 2vrfci_i. (7.13) 



Ci.j Ci^-i 

This corrected equation is precisely the one we would obtain from finite size corrections to 
the su{n) NBA equations. To find this equation from the Bethe ansatz point of view one 
can simply repeat either of the derivations in section |3|, that is the known transfer matrices 
in various representations or the bosonic duality described in the previous sections. In 
section ^ we consider the AdS/CFT Bethe ansatz equations which are based on a large 
rank symmetry group, namely PSU{2,2\A). There one can see an example of how this 
could be done in practice (we will only use the dualities approach because at present we 
don't have the P«S'[/(2, 2|4) transfer matrices for this (exotic) Bethe ansatz equations.). 
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